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ABSTRACT. Three stations near the calving front of the Ross Ice Shelf, Antarctica, recorded GPS data 
through a full spring-neap tidal cycle in November 2005. The data revealed a diurnal horizontal motion 
that varied both along and transverse to the long-term average velocity direction, similar to tidal signals 
observed in other ice shelves and ice streams. Based on its periodicity, it was hypothesized that the 
signal represents a flow response of the Ross Ice Shelf to the diurnal tides of the Ross Sea. To assess the 
influence of the tide on the ice-shelf motion, two hypotheses were developed. The first addressed the 
direct response of the ice shelf to tidal forcing, such as forces due to sea-surface slopes or forces due to 
sub-ice-shelf currents. The second involved the indirect response of ice-shelf flow to the tidal signals 
observed in the ice streams that source the ice shelf. A finite-element model, based on viscous creep 
flow, was developed to test these hypotheses, but succeeded only in falsifying both hypotheses, i.e. 
showing that direct tidal effects produce too small a response, and indirect tidal effects produce a 
response that is not smooth in time. This nullification suggests that a combination of viscous and elastic 
deformation is required to explain the observations. 
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INTRODUCTION 

GPS data collected near the seaward front of the Ross Ice 
Shelf, Antarctica (at stations NN, NS and R13 near a part of 
the Ross Ice Shelf referred to as 'Nascent Iceberg'; Fig. 1), in 
November 2005 indicated that the horizontal ice flow was 
modulated by the ocean tides (Brunt and others, 2010), 
which are strongly diurnal in the Ross Sea (Williams and 
Robinson, 1980; MacAyeal, 1984). This motion, over the 
course of a day, produced a flow rate acceleration and a 
deceleration, or a velocity variation, that was comparable in 
magnitude to the time-averaged annual mean flow (Fig. 2a). 

Similar tidal signals have been observed in other ice 
shelves (Riedel and others, 1999; Doake and others, 2002; 
King and others, 2011; Makinson and others, 2012), ice 
streams (Anandakrishnan and Alley, 1997; Anandakrishnan 
and others, 2003; Bindschadler and others, 2003a, b; 
Gudmundsson, 2006, 2007; Murray and others, 2007), 
outlet glaciers (Marsh and others, 2013) and tidewater 
glaciers (O'Neel and others, 2001). These studies used 
observation techniques to link horizontal motions in the ice 
to a causal mechanism involving the tide based on 
frequency. Most of these studies are either limited with 
respect to geographic extent or restricted to either the ice- 
shelf or ice-stream regime. 

Riedel and others (1999) observed a tide-induced hori- 
zontal displacement of the southern grounding zone of the 
Ekstrom Ice Shelf, in a 1 6 day GPS record (Fig. 1 ). Doake and 
others (2002) used a continuous, 30 day GPS record to 
examine the tide-induced, three-dimensional (3-D) motion 
of the Brunt Ice Shelf. They concluded that the observed 
variations in velocity might have been the result of basal 
friction, or drag, associated with tidal currents interacting 
with the underside of the ice shelf. However, the basal 
friction coefficients needed to justify this conclusion were 
found to be too high to be physically reasonable. 


King and others (2011) collected longer (2 month to 
2 year) GPS records from three sites on the Larsen C ice shelf 
and observed semidiurnal, diurnal and fortnightly modula- 
tions in horizontal ice-flow velocity. They observed that at 
diurnal frequencies, the ice-shelf velocity varied by up to 
±100% from the time-averaged flow. Their model of this 
ocean/ice-shelf system suggested a nonlinear interaction 
between ocean tidal forcing and ice-flow response. 

More recently, on the Ronne Ice Shelf, Makinson and 
others (2012) observed semidiurnal and diurnal velocity 
variations that reached ±300% of the time-averaged flow. 
Their data came from a large campaign that included nine 
sites that sampled at periods ranging from 2 weeks to 1 year. 
Based on tidal phasing, tilt calculations and harmonic 
analysis, they concluded that the velocity variations were 
the result of an elastic response to tide-induced surface 
slope. 

Various studies have observed tidal signals in the 
horizontal flow of the fast-moving West Antarctic ice 
streams that feed the Ross Ice Shelf (Fig. 1). Anandakrishnan 
and others (2003) observed smooth, periodic tidal modula- 
tions in the horizontal flow of Bindschadler Ice Stream (Fig. 
2b). They found that the magnitude of the response decayed 
upstream, away from the grounding zone. Because the 
interaction between the ice shelf and the ice streams was not 
the main focus of their study, Anandakrishnan and others 
(2003) conceded that the influence of the ice shelf on the 
motion of the ice streams was poorly understood. 

Bindschadler and others (2003a, b) and Wiens and others 
(2008) used GPS to investigate a tidal signal in the horizontal 
flow of Whillans Ice Stream (Fig. 2c). However, they 
observed the signal to have a quasi-periodic, 'stick-slip' 
motion. This motion changed from accelerations to decel- 
erations on very short timescales, and was interpreted to be 
associated with subglacial conditions that are sensitive to 
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Fig. 1. Nascent Iceberg GPS stations (NN, NS and R13) on MODIS 
image (Scambos and others, 2007) and InSAR ice velocity (Rignot 
and others, 2011). Black line is MODIS grounding zone. Ekstrom 
(El), Brunt (BISh), Larsen C (LCISh) and Filchner-Ronne (FRISh) ice 
shelves, Rutford (RISt), Mercer (A), Whillans (B), Kamb (C), 
Bindschadler (D) and MacAyeal (E) ice streams, Roosevelt Island 
(Rl), Beardmore Glacier (BG), Cape Hallett (CH) and sites B01 0 and 
DFLT are also indicated. 


the state of the ocean tide in the Ross Sea Embayment. The 
pulsed ice-stream discharge involved a short-term cycle of 
stress accumulation and triggered stress release associated 
with sudden failure of subglacial till, which has been 
modeled as a plastic material (Winberry and others, 2009). 
The GPS measurements indicated that the pulsed ice-stream 
discharge occurred with a duration that is much smaller than 
a day (~20 min). The fact that these pulsations, typically two 
a day, were coordinated in time with the high or falling tide 
in the Ross Sea Embayment suggested that ocean tides 
influenced the dynamics of pulsed ice-stream flow. 

Observations of tidal variations in grounded ice, up- 
stream of the Ross Ice Shelf, are not limited to the fast- 


flowing West Antarctic ice streams. Using GPS data, Marsh 
and others (2013) observed velocity variations on both sides 
of the grounding zone in the vicinity of Beardmore Glacier, 
an outlet glacier which is sourced from the East Antarctic ice 
sheet. Downstream of the grounding zone, the velocity 
varied by 50% from the time-averaged flow. However, once 
grounded, the velocity variation decayed rapidly, within 
15 km of the grounding zone. 

Tidal modulation of ice flow has been observed in 
relatively long GPS records on the Rutford Ice Stream, which 
feeds the Ronne Ice Shelf. Tides in this region are 
predominantly semidiurnal (Padman and others, 2002). In 
a study that extended both upstream and downstream of the 
Rutford Ice Stream grounding zone, Gudmundsson (2006) 
observed diurnal, semidiurnal and fortnightly tidal modula- 
tion in a 7 week GPS survey. Subsequently, Gudmundsson 
(2007) developed a model to investigate the relationship 
between the ice-flow modulation observed on fortnightly 
timescales and the basal shear-stress distribution upstream of 
the grounding zone, and concluded that the strong variation 
on the lunisolar synodic fortnightly tidal constituent (MSf) 
period (14.76 days) was the result of nonlinear interaction 
between the semidiurnal tidal constituents. Further, Gud- 
mundsson (2007) concluded that other ice streams in the 
Filchner-Ronne Ice Shelf system were expected to have 
similar fortnightly flow variations. Ultimately, these obser- 
vations and modeling efforts evolved into a coupled ice- 
stream and ice-shelf nonlinear basal sliding law (Gud- 
mundsson, 2011). In a longer, 750 day GPS survey at a site 
~40 km upstream of the grounding zone, Murray and others 
(2007) reported tidal signals on semidiurnal, diurnal, 
fortnightly, semi-annual and annual periods. In each of the 
Rutford Ice Stream GPS surveys, the observed velocity 
variation was smooth and strongly correlated with the tides, 
similar to the signal observed on Bindschadler Ice Stream. 

In light of investigations on other ice shelves and ice 
streams, the diurnal signal observed in the GPS data 
collected at Nascent Iceberg poses two compelling ques- 
tions: What mechanism, or combination of mechanisms, 
associated with the tide generates the velocity variations 
observed at the front of the Ross Ice Shelf? If there are 
multiple forcing mechanisms, what is their relative signifi- 
cance? 


DATA 

Nascent Iceberg GPS data 

GPS data were collected at the front of the Ross Ice Shelf at a 
site named 'Nascent Iceberg' (78.1 °S, 178.5°W), in close 
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Fig. 2. Time series of ice-flow velocities (solid lines), fluctuating about the time-mean (dashed line), for (a) NN, (b) DFLT and (c) B01 0, based 
on GPS data spanning at least a full spring-neap tide cycle (Brunt and others, 2010). (Data for DFLT and B010, R. Bindschadler, personal 
communication, 2007.) 
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proximity to a large detachment rift, expected to eventually 
yield an iceberg (Brunt and others, 2010). Two GPS stations 
were established, 4.2 km apart, on either side of the rift 
(Nascent North ('NN') and Nascent South ('NS'); Fig. 1). A 
third station (R1 3) was deployed ~50 km inland from the rift. 
At each station, a dual-frequency GPS receiver sampled data 
continuously for at least 1 6 days. 

The GPS data were processed using Precise Point 
Positioning techniques. The residual signal observed at each 
station in the filtered data, which included the removal of 
spurious outliers, and the smoothing and detrending of the 
full dataset, indicated smooth, sinusoidal motions corres- 
ponding to tidal frequencies (Fig. 2a). Furthermore, these 
motions were comparable to, or exceeded, the daily mean 
flow of the ice shelf derived from long-term averages of the 
GPS data (Brunt and others, 2010). 

In comparison to GPS records from the West Antarctic ice 
streams, collected over a ~75 day period during the 2003/04 
field season (personal communication from R. Bindschadler, 
2007), the horizontal variations at Nascent Iceberg (Fig. 2a) 
were similar to the variations observed at Bindschadler Ice 
Stream (Fig. 2b). However, the tide-triggered stick-slip 
motion observed on Whillans Ice Stream was unique, where 
most of the ice flow occurred during two 20 min slip 
intervals over the course of a day (Fig. 2c). 


VISCOUS ICE-FLOW MODEL 

To investigate the possible cause of the tidal variation at 
Nascent Iceberg, a planar, two-dimensional viscous ice-flow 
model was created to examine the interaction between the 
Ross Ice Shelf and ocean tides. Following MacAyeal (1989), 
the physical basis of this model is a simplification of the 3-D 
stress-balance equations, often referred to as Stokes' equa- 
tions, which make use of various assumptions associated 
with the creep of ice, assumed to obey Glen's flow law. 
Because the model employs a power-law creep rheology as 
opposed to an elastic rheology, we refer to it as a 'viscous' 
ice-flow model. The most important assumptions are that 
accelerations are negligible and that horizontal velocities u 
and v are independent of the vertical coordinate, z. These 
assumptions are also known as the shallow-shelf approxima- 
tion. Under these assumptions, the governing equations are 
written in a manner here that includes tidal forcing terms: 
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where x and y are horizontal coordinates, z is the vertical 
coordinate, u and v are horizontal velocities, V is the depth- 
averaged effective viscosity, h is the ice thickness, 77 is a 
perturbation of the free ocean surface, Z is the astronomical 
tide potential, g is the gravitational constant, p\ and p w are 
the densities of ice and sea water, taken to be 917 and 


1028 kg m 3 , respectively, and t x and r y are the basal shear 
stresses acting in the horizontal plane due to ocean currents. 

Terms appearing on the left-hand side of Eqns (1) and (2) 
represent the divergence of the horizontal components of 
the deviatoric stress tensor within the vertical column of the 
ice shelf. All of these terms contain reference to the depth- 
averaged effective viscosity, v, which is defined as a strain- 
rate dependent function that represents Glen's flow law for 
ice (Glen, 1955). Glen's flow law treats ice as a nonlinear 
viscous material and states that ice strain rates are 
proportional to the deviatoric stress raised to the power rt; 
with respect to the creep of ice, n is generally taken to be 3 
(e.g. Glen, 1958; Hooke, 1981). 

The depth-averaged effective viscosity, v, varies with the 
horizontal velocity components, u and v, and with the 
temperature-dependent rate constant, often referred to as 
either the stiffness or flow-law parameter, B. Values for B 
depend on thickness, h, and thus, are spatially variable. For 
this model application, values for B were prescribed based 
on MacAyeal and Thomas (1986). 

Terms appearing on the right-hand side of Eqns (1 ) and (2) 
represent the forces that affect the flow of the ice shelf. The 
first term on the right-hand side of each equation represents 
the effect of the imbalance between sea-water pressure and 
ice pressure, which drives mean, time-independent ice-shelf 
flow. It is common for ice-shelf models to be run with only 
this first term. The other three terms appearing on the right- 
hand side of Eqns (1 ) and (2) pertain to time-dependent tidal 
forcing. In the first case, sea-surface slopes are induced by 
the deformation of the free surface of the ocean, r,, as the 
tidal perturbation passes beneath the ice shelf. In the second 
case, the astronomical tide-generating potential, Z, gener- 
ates a horizontal force that pulls on the ice ever so slightly 
(Thomas, 2007). The last term on the right-hand side of each 
equation pertains to basal shear stress, r, associated with 
tide-induced ocean currents within the ice-shelf cavity. 

Based on these governing equations, a finite-element 
model of the Ross Ice Shelf was created. The numerical 
domain of this model was a mesh composed of 1500 nodes 
that create 2700 triangular elements. The domain boundary 
of the model was based on the NASA Moderate Resolution 
Imaging Spectroradiometer (MODIS) Mosaic of Antarctica 
(Scambos and others, 2007). Where the domain boundary 
was in direct contact with grounded ice (e.g. ice streams or 
outlet glaciers), the associated nodes were designated as 
regions of 'inflow', and the horizontal velocities at these 
nodes were specified based on interferometric synthetic 
aperture radar (InSAR) methods of Rignot and others (201 1) 
(Fig. 3a). This designation accounted for 20% of the 
boundary. Where the domain boundary was in direct 
contact with rock (e.g. adjacent to mountains), the hori- 
zontal velocity of the ice shelf at these nodes was initially 
specified to be zero, and conceptually represents 'no flow' 
across the grounding zone; this accounted for 60% of the 
boundary. Finally, on the seaward front of the domain, the 
total stress acting in the horizontal plane at the ice-front cliff 
was specified to be the vertically integrated pressure in the 
sea water beyond. These nodes were designated as 'outflow' 
and accounted for 20% of the boundary. 

Ice thickness, h, varies spatially within the Ross Ice Shelf 
and was estimated based on the densities of ice and water, 
using British Antarctic Survey BEDMAP surface elevation 
data (Lythe and others, 2001) and assuming hydrostatic 
equilibrium. A constant value of 20m was removed to 
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Fig. 3. (a) Velocity based on InSAR (Rignot and others, 201 1) compared with (b) model-derived velocity, with finite-element mesh (black 
triangles) for reference. 


approximate a firn correction. Ultimately, thickness results 
based on hydrostatic assumptions compared well to radio- 
echo sounding estimates of Neal (1979). 

The finite-element model was then implemented in three 
distinct patterns of usage. In the first pattern, the steady-state, 
time-independent, ice-shelf flow was solved for under the 
assumption that time-dependent tidal forcing (e.g. the last 
three terms on the right-hand side of Eqns (1) and (2)) was 
absent. This model result was then used as a baseline to 
isolate the effects of the model results of the other two 
patterns of usage, which employ either direct or indirect 
tidal forcing. In the second pattern of model usage, the 
'direct tidal forcing model', Z was assumed zero (because of 
its small magnitude) and 77, t x and r y in Eqns (1) and (2) were 
extracted from an ocean tidal model (Padman and others, 
2003) and specified at each node through a full spring-neap 
tide cycle to produce a solution that varied through time. In 
the third pattern of model usage, the 'indirect tidal forcing 
model', the terms Z, 77, r x and r y were assumed zero, but the 
ice flow along 'inflow' boundaries was specified to have a 
tidal variation, or stick-slip variation, depending on the 
model experiment. This means that the difference between 
the second and third patterns of model forcing was designed 
to differentiate the effects of forcing by the ocean acting 
directly on the ice shelf from the effects of the ice streams 
acting indirectly on the ice shelf. Ultimately, the Nascent 
Iceberg GPS data were used to interpret the model results of 
all three patterns of usage. 

Steady-state ice-shelf flow 

The first pattern of model usage, defined as the steady-state, 
time-independent, ice-shelf flow, was solved for under the 
assumption that time-dependent forcings were absent. Thus, 
the three terms appearing on the right-hand side of Eqns (1) 
and (2), that pertain to time-varying tidal forcing, were taken 
to be zero. 

To ensure satisfaction of the power-law constitutive 
relation, the model iteratively solved for viscosity and 
horizontal velocity over the entire ice shelf using a 
prescribed initial value for depth-averaged effective viscos- 
ity, V, and prescribed initial horizontal velocity components, 
u and v, in regions of 'inflow' based on Rignot and others 
(2011). With an initial depth-averaged effective viscosity of 
17=1 x10 17 Nsm~ 2 , 25 iterations were needed for the 
model to converge on solutions for viscosity and velocity. 


The model results for v ranged from 10 13 to 10 17 Nsrrr 2 , 
with a mean value of 9.8 x 1 0 14 N s rrr 2 . 

The model result for steady-state, time-independent, ice- 
shelf flow (Fig. 3b) confirmed that equations, input data, 
mesh resolution and coding were consistent with the mean 
state of ice-shelf flow. Specifically, model results for velocity 
were compared to the InSAR velocities of Rignot and others 
(201 1) (Fig. 3a) and the time-mean ice-shelf velocity based 
on GPS data collected at Nascent Iceberg. The model was 
tuned solely using the 20 m approximated firn correction, 
and ultimately derived a horizontal velocity at Nascent 
North of 1 093.87 ma 4 , comparing well to the value of 
1 093.59 ma 4 based on the GPS data. 

Direct, tide-forced, ice-shelf flow model 

Tides of the Ross Ice Shelf are strongly diurnal (Williams and 
Robinson, 1980; MacAyeal, 1984), with maximum tidal 
amplitudes occurring near the grounding zones, specifically 
in the vicinity of Mercer Ice Stream (Fig. 1), where the total 
amplitude range is >3 m (Padman and others, 2003). The 
tidal perturbation rotates clockwise under the floating Ross 
Ice Shelf like a Kelvin wave, around a diurnal amphidromic 
point located off the coast near Cape Hallett (MacAyeal, 
1984) and a weaker semidiurnal amphidromic point 
beneath the northwestern part of the Ross Ice Shelf (Williams 
and Robinson, 1980). 

The tidal amplitude, or perturbation of the free surface, 77, 
in Eqns (1) and (2), represents small changes in height across 
the ice shelf. Given the scale of the Ross Ice Shelf (~800 km 
on each side), at any one time, some regions are experien- 
cing high tide, while others are experiencing low tide. 
Variable heights across the ice shelf will create regional 
surface slopes; these slopes may be substantial enough to 
induce accelerated and decelerated ice flow. The term in 
Eqns (1) and (2) containing Z describes the astronomical 
tide-generating potential, the effect of which is equivalent to 
putting the ice shelf on an inclined ocean surface, induced 
by the equilibrium tide. Based on Platzman (1984), the 
equilibrium tide has an amplitude of 0.1 m and, due to the 
diurnal nature of the tides in the Ross Sea Embayment, a very 
large phase variation, equivalent to half the distance around 
the circumference of the globe, which is of order 1 0 000 km 
at 75-85° S. In contrast, sea-surface tidal amplitudes for the 
Ross Sea Embayment are on the order of 1 m and vary on 
length scales of 1000 km. Thus, the theoretical slope 
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Fig. 4. An element-by-element difference between the velocity baseline (Fig. 3b) and (a) the influence of r, (b) r\ and (c) the combined 
influence of r and r\, representing the direct, tide-forced velocity. 


associated with the term containing Z was disregarded as it 
is ~100 times smaller than the term associated with sea- 
surface slope, containing 77 . The last term in Eqns (1) and (2), 
containing r, is associated with basal friction resulting from 
the interaction between the base of the ice shelf and the 
ocean currents within the ice-shelf cavity. Robertson and 
others (1998) define the basal shear stresses due to ocean 
currents in the x and y directions within the ice-shelf cavity, 
t x and T y , as 

Cd Pw U w | U w| \ 

Tx = Tp (3) 

_ Co Pw V w \U 

T ' = H2 ' ( 4 ) 

where Co is the coefficient of friction taken to be 0.003 from 
Robertson and others (1998), p w is the density of water, U w 
and V w are the horizontal tidal transports within the ice-shelf 
cavity, U w is the magnitude of the volume water transport 
vector ( | Uw | = v/Ul + V7J), ar| d H is the water column 
thickness under the ice shelf. 

Equations (1 ) and (2) required the specification of 77 and r, 
through time, at each node in the finite-element mesh. The 
free surface, 77, was extracted directly from the Ross Sea 
Height-Based Tidal Inverse Model ('Ross_lnv_2002') (Pad- 
man and others, 2003), which represents a high-resolution 
finite-difference model for the Ross Sea, ideal for predicting 
heights for the Ross Ice Shelf. The ocean tide model has 
~ 10 km resolution and its accuracy is optimized by the 
assimilation of data from ten gravimetric tide gauges on the 
Ross Ice Shelf. Accuracy is likely to be at the cm level for 77 . 
Similarly, to determine the basal shear stress, r, the water 
column thickness, H, which is based on the seismic and 
radar analysis, and horizontal tidal transports within the ice- 
shelf cavity, Uw and V w , were also extracted from the ocean 
tide model and then applied to Eqns (3) and (4). 

Once variables from the tide model were extracted, the 
direct tide-forced, time-dependent, ice-shelf flow model was 
iterated to solve for viscosity and velocity. This process was 
repeated, with a time step of 2 hours, for the 1 7 day length of 
the GPS record, or through slightly more than one spring- 
neap tidal cycle. 

The influence of tide-induced basal friction, or drag (r), 
and that of tide-induced sea-surface slope (associated with 
77 ) on the flow of the ice shelf were initially examined 
separately. To uniquely examine the effect of drag on the 
flow of the ice shelf in the absence of forcing associated with 
slope, a series of model velocity fields was generated for the 


entire Ross Ice Shelf, with one field associated with each 
2 hour time step. A 'maximum drag-forced model velocity 
field' was created by extracting the maximum drag-forced 
model velocity of each element, independent of when it 
occurred during the tidal cycle. The steady-state solution 
(Fig. 3b) of each element was then subtracted from this 
'maximum drag-forced model velocity field' to isolate the 
maximum velocity fluctuation associated with r (Fig. 4a). 
The effect of tide-induced sea-surface slope (associated with 
77 ) on the flow of the ice shelf, in the absence of forcing 
associated with drag (r), was similarly isolated (Fig. 4b). 

A series of model velocity fields to test the collective 
influence of both r and 77 was generated to ultimately assess 
the direct, tide-forced ice-shelf flow model. A 'maximum 
model velocity field' was created by extracting the max- 
imum model velocity of each element, independent of when 
it occurred during the tidal cycle. The steady-state solution 
(Fig. 3b) of each element was then subtracted from this 
'maximum model velocity field' to isolate the maximum 
velocity fluctuation associated with both tide-induced sur- 
face slope and basal drag (Fig. 4c). 

From Figure 4a and b, the influence of r, which peaked at 
2.0 m a -1 , was an order of magnitude less than the influence 
of 77 , which peaked at 41.4ma~ 1 . Further, the forcing 
associated with 77 is symmetric, while the influence of r was 
not. The maximum response of both sea-surface slope and 
basal drag occurred along the front of the ice shelf. For the 
slope forcing, this occurred in the vicinity of Nascent 
Iceberg. However, the maximum response associated with 
drag forcing moved along the front of the ice shelf 
depending on time. After high tide, the maximum response 
occurred in the vicinity of Roosevelt Island and then 
migrated toward Nascent Iceberg. This migration was likely 
due to the fact that, as discussed previously, the tidal 
perturbation rotates clockwise under the floating Ross Ice 
Shelf around an amphidromic point. 

Based on model results, the collective effect of both tide- 
induced surface slope and basal drag (Fig. 4c) was not 
enough to account for the variations observed at Nascent 
Iceberg. The maximum contribution from the direct tidal 
forcing was taken to be the difference between ( 1 ) the 
maximum value of each element over the full spring-neap 
tidal cycle and (2) the baseline velocity (Fig. 4). The 
maximum difference, which occurred in the vicinity of 
Nascent Iceberg, was ~42.9ma~ 1 , or <5% of the 1000m 
a -1 velocity variation observed in the GPS data (Fig. 2a). 
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Fig. 5. (a) (7(f) observed at DFLT (red line), where U (black solid line) and ma x((7(t)) (black dashed line) are used to prescribe a sine wave to 
the ice-stream forced flow model, (b) A detailed view of the period 19-21 lanuary 2004 (shaded gray box in (a)) indicates the form of the 
prescribed sine wave (blue line). 


Thus, based on this viscous flow model, tide-induced 
surface slope and basal drag were insufficient to explain 
the magnitude of the velocity fluctuations observed at the 
front of the Ross Ice Shelf. 

Indirect, ice-stream forced, ice-shelf flow model 

To examine the indirect tidal influence on the ice-shelf flow 
model, i.e. forcing associated with a time-dependent vari- 
ation in 'inflow' from the ice streams, Eqns (1) and (2) were 
solved for under the assumption that time-dependent 
forcings were absent with a modified set of boundary 
conditions. Specifically, the velocities of regions of 'inflow' 
associated with the ice streams were allowed to vary through 
time based on the magnitude and character (either 
sinusoidal or stick-slip) of the West Antarctic ice-stream 
GPS data (Fig. 2b and c). The velocities of regions of 'inflow' 
associated with the outlet glaciers were specified based on 
the InSAR velocities of Rignot and others (2011) and held 
constant. 

The full ice-flow velocity, (7(f), unique to each ice 
stream, is defined as 

17(f) = U+ U'(t) , (5) 

where U is the time-independent, mean velocity of the ice 
stream, and (/(f) is the time-dependent part of the velocity 
field that is dominated by the tide (Fig. 5a). The mean 
component, 17, was calculated from the ice-stream GPS data 
using the least-squares method. The magnitude of U in the 
vicinity of the grounding zone varied between ice streams 
from near-zero to ~400 m a~’ . The time-dependent velocity, 
I7'(f), was tabulated using 21 tidal constituents based on a 
similar study on Rutford Ice Stream by Murray and others 
(2007). Given that the ice-stream GPS dataset was collected 
over ~7 5 days, long-period tidal variation is evident in 
I7'(f). The magnitude of I7'(f) varied substantially between 
the ice streams from near zero to ~1 0 000 m a -1 , depending 
on the character of the ice-stream flow (either sinusoidal or 
stick-slip; Fig. 2b and c, respectively). 

Bindschadler Ice Stream was observed to have a smooth, 
diurnal, sinusoidal 17' (f) (Fig. 2b). However, observed U at 


its grounding zone was extremely high, resulting in a large 
17(f). The character of the velocity variation observed in this 
ice stream matched the diurnal, sinusoidal motion observed 
at Nascent Iceberg; however, the magnitude of the variation 
was amplified at the ice-shelf front (Fig. 2a and b). Whillans 
Ice Stream was observed to have a tide-triggered, stick-slip 
I7'(f). The slip events were very brief (two 20 min slip 
intervals per day), but during those slip events, velocity 
observations were on the order of lOOOOma -1 ; velocity 
between slip events was observed to be approximately zero 
(Fig. 2c). 

For this model experiment, the 'inflow' velocities for 
Mercer, Karnb, Bindschadler and MacAyeal ice streams were 
varied sinusoidally based on the ice-stream GPS data. 
Specifically, the mean velocities, U, were modified by 
defining a simplified sinusoid around that value (Fig. 5b), 
with a maximum amplitude limited by the percent 
change between 17(f) and the maximum full velocity 
derived from the GPS data, unique to each ice stream, 
i.e. max( 17(f)) = U+ max(t7'(f)). Observations of tide- 
triggered, stick-slip 17' (f) have thus far been limited to 
Whillans Ice Stream (Bindschadler and others, 2003a, b). 
Therefore, stick-slip motion in the model was only 
prescribed to this ice stream. Values of 'inflow' associated 
with the outlet glaciers were based on the initial InSAR 
velocities and held constant. 

Based on model results, the indirect ice-stream forcing 
was not enough to account for the variations observed at 
Nascent Iceberg. The maximum contribution from the 
indirect ice-stream forcing was taken to be the difference 
between (1) maximum model velocity of each element over 
all 'time steps' and (2) the baseline velocity (Fig. 6), and is 
greatest in the vicinity of Whillans Ice Stream. At the front of 
Ross Ice Shelf, this value was 580 ma -1 , or ~50% velocity 
fluctuation observed in the GPS data collected at Nascent 
Iceberg (Fig. 2a). Relative to the results of the surface slope 
and basal friction model, the magnitude of the response of 
the ice-stream forced model was considerably larger. 
However, the model velocity field at Nascent Iceberg spiked 
in response to the stick-slip model forcing at the Whillans 
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Fig. 6. The difference between the indirect, stick-slip (for Whillans 
Ice Stream only) and sinusoidal (all other ice streams) stream-forced 
model velocity maximum of each element and the velocity baseline 
(Fig. 3b). 


Ice Stream region of 'inflow'. Thus, the character of the 
model response at Nascent Iceberg did not match the 
smooth, periodic observations of the GPS record. 



1060 I I I I I I 

Diurnal tide cycles 


b 



Fig. 7. The magnitude and character of model velocity at Nascent 
Iceberg for (a) direct tidal forcing (response to surface slope and 
basal friction) and (b) indirect tidal forcing (response to changes in 
'inflow' from the ice streams) for at least one tidal cycle. Note the 
changes in the y-axes. 


DISCUSSION 

The two sets of experiments conducted here have failed to 
simultaneously predict the magnitude and the smooth, 
periodic nature of the observed variations in ice velocity at 
the front of the Ross Ice Shelf. The model results associated 
with tide-induced surface slopes and basal friction generally 
reproduced the smooth, diurnal nature of the observations at 
Nascent North (cf. Figs 2a and 7a); however, this model 
application does not reproduce the magnitude of the 
observations. Conversely, the model results associated with 
tidal variation of ice-stream 'inflow' velocities generally 
produced larger-magnitude variations, relative to the direct, 
tide-forced, ice-shelf flow model (cf. Figs 2a and 7b); 
however, this model application does not reproduce the 
smooth character of the observations. To elicit the smooth, 
sinusoidal variation of the ice-shelf flow at Nascent Iceberg, 
the discharge across the grounding line of Whillans Ice 
Stream would also have to be smooth and sinusoidally 
varying, not stick-slip. Also, the amplitude of this discharge 
would have to be much higher than observed, and higher 
even than the amplitude at Nascent Iceberg, because the 
amplitude pattern decays downstream of the input bound- 
aries. Ultimately, in the case of tidal forcing due to sea- 
surface slopes or due to sub-ice-shelf currents, the viscous 
stress regime of the ice shelf is too stiff to allow an 
appropriately large strain rate demanded by the GPS 
observations at Nascent Iceberg. These results nullify the 
two hypotheses posed above, but yield a firm conclusion 
that 'viscous'-style flow responses to the tidal forcing, either 
indirect or direct, must be augmented by other rheological 
processes to explain the GPS data. 

Harmonic analysis of GPS data from the Ronne Ice Shelf 
(Makinson and others, 2012) suggested that tide-induced 
surface slope alone was capable of producing the observed 
±300% semidiurnal and diurnal velocity variations around 
the time-averaged ice-shelf flow. Their conclusions also 


minimized the effect of basal drag, mostly based on phasing 
arguments. Further, based on the temperature of the ice and 
modeling efforts of Gudmundsson (2007), they concluded 
that the Ronne Ice Shelf responds as an elastic membrane at 
tidal timescales. 

Given the results of Makinson and others (2012), and 
given that the character of the model results associated with 
tide-induced surface slopes and basal friction more closely 
matched the smooth nature of the velocity variations at 
Nascent Iceberg, we conclude that this model provides the 
viscous component of a viscoelastic ice-flow model, which 
assumes that ice behaves as a Maxwell material. The 
smooth, periodic nature of the signal observed at Nascent 
Iceberg further suggests that some component of elastic 
response must be dampening the stick-slip signal observed 
at Whillans Ice Stream. 

The timescale, f, of viscoelastic relaxation of a Maxwell 
material is characterized by 


where E is the Young's modulus, taken to be 4.8 GPa based 
on modeling efforts of Gudmundsson (2011), and v is the 
effective viscosity, taken to be 9.8 x 10 14 N srrr 2 based on 
the mean depth-averaged effective viscosity, v, determined 
by the Ross Ice Shelf flow model. The relaxation time is 
therefore on the order of a day, or roughly tidal timescales. 
Therefore, we conclude that the inclusion of a viscous 
component, and the incorporation of an elastic component, 
in our model is appropriate for assessing ice-flow response 
to tidal forcing. 

For a Maxwell material, the total strain of the system 
(etotai) is the sum of the viscous and elastic strains (eviscous and 
^elastic) on eac h element: 

Etotal Eviscous ± ^elastic. ( w 

Under this assumption, the viscous ice-flow model pre- 
sented here has determined the component associated with 
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the viscous response (e V iscous); the elastic component Elastic) 
can now be solved independently, or in series, with respect 
to the viscous component. Ultimately, we conclude that 
adding the elastic component to the existing ice-shelf flow 
model will capture the magnitude of the signal observed at 
Nascent Iceberg while simultaneously explaining the down- 
stream loss, or dampening, of the signal observed at 
Whillans Ice Stream. 

CONCLUSIONS 

Two weeks of CPS data were collected at the front of the Ross 
Ice Shelf and contained a diurnal signal, assumed to be 
associated with the tide. A viscous ice-flow model, using the 
nonlinear viscosity implied by Glen's flow law, was created 
to assess the relative significance of direct and indirect tidal 
forcing mechanisms based on the response observed at the 
front of the Ross Ice Shelf. Specifically, model experiments 
were designed to independently quantify the effects of (1) 
tide-induced surface slopes and basal friction, and (2) tide- 
modulated flow of the ice streams. In general, the model 
response to tide-induced surface slopes and basal friction 
matched the character of the observations at Nascent 
Iceberg, but was too small in amplitude. While the response 
to tide-modulated flow of the ice streams more closely 
matched the magnitude of the observations at Nascent 
Iceberg, the temporal variation of the response in the model 
was not smooth enough to match the observations. These two 
experiments suggest that there is no possible way to force a 
viscous treatment of the ice shelf to match observations, and 
that some additional aspect of ice rheology must come into 
play in order to explain the observations. 

Ultimately, the effect of the tide-induced slope of the ice 
shelf, coupled with the effect of basal drag associated with 
tidal currents in the ice-shelf cavity, amounted to <5% of the 
horizontal velocity variation observed at Nascent Iceberg. 
However, based on the analysis and conclusions of 
Makinson and others (2012), this model represents the 
viscous component of a viscoelastic ice-flow model, where 
the solution for the elastic component can be added to this 
result, under the assumption that ice behaves as a Maxwell 
material. Modeling and understanding the periodicity and 
sensitivities of ice-shelf flow will illuminate factors associ- 
ated with the long-term stability of ice shelves. Further, a 
better understanding of the mechanisms controlling the tidal 
signals on this ice shelf may provide insight into the 
mechanisms at play across the grounding zone, within the 
complex regime of the ice streams. 
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